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C^- I Abstract 

^vq , We present a multicanonical algorithm for the SU(3) pure gauge theory at the deconfinement 

phase transition. We measure the tunneling times for lattices of size L 3 x 2 for L = 8, 10, and 
12. In contrast to the canonical algorithm the tunneling time increases only moderately with 
L. Finally, we determine the interfacial free energy applying the multicanonical algorithm. 

£2 ■ 1 Introduction 

For the analysis of systems with first order phase transitions-as e.g. quenched QCD at the 
deconfinement temperature-surface effects are of crucial importance. They determine the 
physical dynamics and lead to supercritical slowing down in computer simulations at the 
critical point. Configurations which include an interface are exponentially suppressed by 
their free energy which results in an exponential increase of the tunneling time. Recently 
the multicanonical algorithm was proposed in order to overcome this supercritical slowing 
down [|l]. Here we present a variant of this method for the SU(3) gauge theory with the 
Wilson action. We have used lattices of size L 3 x Lt with Lt = 2 and L = 8, 10, and 12 for 
which there is a strong first order phase transition. We achieve a considerable reduction in 
the tunneling time with only a modest increase in computer time. We apply this algorithm 
to the determination of the interfacial free energy 

F cd = a cd A (1.1) 

of a confmcment-dcconfinemcnt interface of area A, a c d being its surface tension. For this 
we made use of Binder's histogram method which in contrast to the methods used in ^ || 
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generates interfaces dynamically rather than by introducing artificial external parameters. 
In that respect, this approach is similar to transfer matrix methods Q. 

2 Monte Carlo Simulations at a First Order Phase Tran- 
sition 

As mentioned above, MC simulations close to first order phase transitions suffer from super- 
critical slowing down, i.e. the tunneling time tx for a system of volume V — L 3 is expected 
to be proportional to the inverse of the probability of a system with interfaces (actually 
there will be two interfaces because of the boundary conditions) , 

t l cxexp(2L 2 |) (2.1) 

and therefore diverges exponentially with the area A = L 2 of an interface. 

In order to overcome this problem, the multicanonical algorithm does not sample the 
configurations with the canonical Boltzmann weight 

P£ a ™(S)cxexp(/3S), (2.2) 

where S = 1/3 Y] n tr Ua is the Wilson action in four dimensions, but rather with a modified 
weight 

Vr(S) oc eMPL(S)S + a L (S)). (2.3) 

The coefficients and (3l are chosen such that the probability Pl (not to be confused with 
the Boltzmann weights) is increased for all values of the action in between the two maxima 
S™ 1 ' 1 and 5 , ™ aa: ' 2 } as shown schematically in figure 1. Finally data are analyzed from the 
multicanonical samples by reweighting with exp((/3 — Pl{S))S — az,(5)). 

In order to approximate the weights 'P£ lc leading to the distribution P™ c of figure 1 , we 
start from some good estimate of the canonical distribution Pjr (see below) at the coupling 
/3£ which corresponds to equal weight in both phases. This will lead to different heights of 
the two peaks. Then we take a partition 

S° L = -6L 3 L t < Si = S™ axA <Sl<...<S* /2 = Sr n <■■ 

... < = S™ ax ' 2 < — 6L 3 Lt (2.4) 

of the interval —6L 3 L t < S < 6L 3 L t . The coefficients and (3^ are chosen to be constants 
a k L and fi\ in the intervals \S\, S^ +1 ) such that lnP™ c interpolates the linear function 

o nmax.l 

ln praax,l ( m pmax,2 _ fopnax.U . g ~ ,3 5) 
L V Xj n 1 ^max,2 gmax,l v / 

continuously between the points and S k+1 for k = 1 to N — 2. The Boltzmann weight is 
identical to the canonical one in the first and last interval. We arrive at 

,k = 

01-0=1 S0 + ln( 1 ^^)/(S k L +1 -S k L ) ,k = l,...,N-2 (2.6) 

,k = N — 1 
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where Sf3 = In 



max.2 



gmax,i^ an( j ^ e Q,fe are gi ven by 



, k+l 



6p-S 



max,l 
L 



(2.7) 



The partition {S*} which is in principle arbitrary is defined by demanding 



Pl(S : 



1/n ,ft=l,...,JV/2- 
r 2 , k — N/2, . . . ,N 



(2. 
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max, pmax,2 ^ introducing a (5/3 ^ 0. In addition the specific 



ymax,2 
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rmn 
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This generalizes the formulas 

given in to the case P- 
choice for r\ and assures that the probability of S™ in is lifted by the correct amount. 

We apply this algorithm to the SU(3) pure gauge theory at the deconfinement phase 
transition. The multicanonical data sampling was done with a 5-hit Metropolis algorithm 



as well as with a Creutz heat bath algorithm modified according to equ. (2.3). In both cases 
three independent S'C/(2)-subgroups are updated following the idea of ref. ||. Because of 
the dependence of and 0l on the total action S, the update of the active link has to be 
done in scalar mode. However, most of the computation time is needed for the calculation 
of the staples surrounding the active link, and this is still vectorizable. In addition, for 
the Creutz algorithm we profited from the fact that as long as the action stays within one 
interval [S k ,S^ +1 ) there is no dependence of oll and 0l on S. Therefore, if the intervals 
are not too small even the update can be vectorized partially, since one hardly ever crosses 
an interval boundary. 

In table || the update times for the various algorithms for one processor of a CRAY- 
YMP are given for a 10 3 x 2 lattice using N = 14 intervals for "P™ c . As can be seen the 
CPU time per link and update is increased by only about 50%. In our simulations we use 



Table 1. CPU time per link update for canonical and multicanonical 
versions of a Metropolis and a Creutz heatbath algorithm. 



Algorithmus 


Metr. 


Metr.(mc) 


Creutz 


Creutz(mc) 


t U pd{psec) 


37 


53 


31 


52 



C-periodic boundary conditions in the three spatial directions (see ||). This suppresses 
the two deconfined phases with nonvanishing expectation value of the imaginary part of 
the Polyakov line. Then only the confined and one deconfined phase coexist at the critical 
temperature T c . For finite volumes we define T£ (or /3£ resp.) by demanding equal weight 
for these two phases [Q, i.e. 

P L (S < S r £ m ) = P L (S > Sr n ). (2.8) 

We typically do 50000 updates for One sweep consists of one Metropolis (resp. 

Creutz) update step plus four overrelaxation steps. For the smaller lattices L — 8 and 10 
the parameters for the multicanonical algorithm were taken from canonical simulations. For 
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L = 12 where the canonical algorithm is no longer efficient due to the increasing tunneling 
time, we took the parameters from an appropriate finite size scaling extrapolation of the 
distributions for the 10 3 x 2 lattice. In all cases we take N = 14 intervals. 

Figure 2 shows part of the history of the average Polyakov line as a function of the number 
of sweeps for the various algorithms. From the full histories we extract the tunneling times 
7X, defined as one fourth of the number of sweeps needed for a full tunneling cycle. The 
tunneling times tl in table || indicate that the increase with L for the standard algorithm 
is effectively reduced. Already for L = 10 the gain in 17, is much larger than the overhead 
in CPU time per update which is due to scalar operations. 

Table 2. Tunneling times T£ for the various algorithms and lattices of size i 3 X 2. 



L 


a 


Mctr. 


Metr.(mc) 


Creutz(mc) 


8 


5.094 


300 


200 




10 


5.095 


1500 


400 


400 


12 


5.0928 




700 


700 



lattices of 

Thus the multicanonical algorithm provides an efficient tool for Monte Carlo studies of 
the deconfinement phase transition. As an application, we present in the next section the 
determination of the interfacial free energy Fg ■ 



3 The Interfacial Free Energy 

The configurations with S = S™ m are expected to contain two parallel interfaces of area 
L 2 between the confined and the deconfined phase. Therefore they should be suppressed by 
the interfacial free energy aL 2 . We expect additional power law corrections due to capillary 
wave fluctuations and translational Goldstone modes of the interfaces, such that 

Pr"°cL*exp(-2L 2 |) (3-1) 

(see also ref.H). The probability distribution should be flat around S™ m for non- interacting 
surfaces due to translations of the interfaces relative to each other. 

We extract the probability distribution Pj, for the plaquette and the real part of the 
Polyakov line from our multicanonical simulations for L = 8, 10, and 12 (see figure 3). They 
agree reasonably with those obtained in ref. for periodic boundary conditions. However, 
there is no plateau for the probabilities around S™ in . This is likely to be due to interfacial 
interactions which will be reduced by using a larger extension in one of the spatial directions. 
This is left for further investigations. 

Nevertheless, we consider the two quantities 
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F £) = _L Ih^-t- (3.2) 

2L prmn V / 
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and 

Ff^-^lnP™" (3.3) 

where p L max = ifpmax, ^ pmax, y Xhey should converge to y for L — > oo for non- 
interacting surfaces. We plot them as functions of 1/L 2 in figure 4. We find a reasonable 
agreement between the values taken from the plaquette and Polyakov line distributions. As 
expected, the discrepancy between F L ^ and decreases for increasing L. Still, compared 
to the values of a obtained by different methods ^ [| the extrapolated result of our data is 
too small. Simulations on asymmetric lattices should clarify this point. 

In conclusion we have demonstrated the power of the multicanonical algorithm for SU(3) 
pure gauge theory. The tunneling time t£ is reduced considerably while the update time is 
only increased by less than 50% compared to standard Metropolis or heat bath algorithms. 
This reduces the computer time needed for investigations of the deconfinement phase tran- 
sition like the one presented here for the interfacial free energy. Using cylindrical lattices 
will probably reduce the interactions between the interfaces thereby allowing the extraction 
of the interface tension. 
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4 Figure Captions 



Figure 1 Schematic representation of of the probability distributions Pj, (straight line) and 
P™ c (dashed line) corresponding to the canonical Boltzmann weight "P£ an and the 
multicanonical one V n L lc . For S < S™ ax ' L and S > S™ ax <, P L and P^ c (normalized 
to their maximum value) are identical. 

Figure 2 Evolution of the average Polyakov line as a function of the total number of sweeps 
for a 10 3 x 2 lattice for (A) a canonical Metropolis algorithm, (B) a multicanonical 
Metropolis algorithm, and (C) a multicanonical Creutz algorithm. 

Figure 3 Distribution for the average plaquette and the real part of the Polyakov line for lattice 
sizes L = 8, 10 and 12. from the average versus InL/L 2 . indicated by a diamond, 

Figure 4 The free energies and from the average plaquette and the real part of the 
Polyakov line versus 1/L 2 . The value for a/T from reference Q is indicated by a 
diamond, the value from reference Q by a filled square, and the value from reference 
by a triangle. 
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